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The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Numerical Investigation of the 2D JV — 2 Wess-Zumino Model 



Christian Wozar 



1. Introduction 

The two-dimensional JV = 2 Wess-Zumino model in the continuum shows no spontaneous su- 
persymmetry breaking. However, a lattice formulation must break (part of) the supersymmetry 
explicitly due to the failure of the Leibniz rule on the lattice. In this model high statistics on 
large lattices are available and supersymmetry restoration effects can be analyzed numerically. The 
restoration of supersymmetry in the continuum limit must be treated carefully as has been demon- 
strated in supersymmetric quantum mechanics [jl|, ^]. One possible way to circumvent relevant 
supersymmetry breaking operators is the application of a blocking transformation to a free theory 
[Q] leading to solutions similar to the Ginsparg -Wilson relation for the chiral symmetry 

A further suggestion keeps at least one supersymmetry on the lattice preserved and goes under 
the name of "Nicolai improvement" []5p. In former works ^ [7|] such improved models using 
Wilson fermions were simulated, and discrepancies to the perturbative result as well as problems 
with the extraction of masses occurred at stronger couplings. In this work (see [Q] for detailed 
analyses) the effects of the Nicolai improvement in the intermediate coupling regime are analyzed 
and compared to results of unimproved simulations. Additionally different fermion formulations 
(standard/twisted Wilson, SLAC) are explored. 

2. The model 



The continuum action with complex field 
<p = (pi + i(p2, 



d 2 x (2dyd(p + ±\W'((p)\ 2 + \j/M\i/) 
M = fd + fd + W"P + + W"P- 

(2.1) 

is invariant under four real supercharges. 
Taken together they satisfy the jV = (2,2) 
superalgebra, and it has been argued that 
one supersymmetry can be preserved on 
the lattice [0]. 

The holomorphic superpotential here 
(see Fig. |) 
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\m(p 2 + ±g(p 3 



(2.2) 



Figure 1: Classical potential |W'(<pi)| from (J2JJ) for van- 
ishing imaginary part (q>2 = 0). In the free theory limit 
(g — > 0) the left minimum is pushed towards minus infinity. 



contains a mass parameter m and defines a dimensionless coupling X = £ This theory possesses a 
discrete 2,\ symmetry which in general is partially broken by a lattice discretization. A perturbative 
expansion in orders of A around the free theory at A = is possible and is used below. 



2.1 Nicolai improvement 

One real supersymmetry can be preserved on the lattice by using the action 



(2.3) 
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in terms of the Nicolai variable £, x = 2(d(p) x + W x with W x = W'(<p x ), Wxy := dW x /d<p y and 

M„=( W ? 2 lA = (^^\. (2.4) 

In terms of the original field <p, the lattice action reads 

5 = £ (2 (30) Xd(p) x + \ \W x f + W x (d(p) x + W x (d(p) x ) +Y,%M xy y r (2.5) 



x xy 

This action only differs from a straightforward discretization by discretized surface terms 

AS = ^(w x (d(p) x + W x (d(p) x ^ (2.6) 

which must vanish in the continuum limit. 
2.2 The lattice discretization 

To preserve the full supersymmetry of the free theory the same lattice derivatives for bosonic and 
fermionic degrees of freedom must be used. In this extensive study we compare three different 
lattice derivatives: 

• Symmetric derivative (djl) = \ (5*+/^ — ?>x-ix,y) with standard Wilson term W x = W'((p x ) — 
I (A<p) x using r = 1 . The Wilson term must be added to W x (and not to the derivative) in order 
to obtain an antisymmetric matrix (d^)^,. This results in a fermion matrix 



1 2S„ w"(if,)S„) 2 



• Symmetric derivative d s with twisted (imaginary) Wilson term W x = W'((p x ) + j(A(p) x re- 
sulting in 

_(w»{< Px )8 xy 24. \ r 

Mxy -{ idxy w"m8 xy ) +y ^- (2 - 8) 

The choice r = 2/\/3 reproduces the mass of the free theory up to G(a A ) at lattice spacing a 
as discussed in 

• SLAC derivative d x ^ y = (— 1 ) x ~ y ^$~yjjm » ^ = ^ w ^ t ' 1 f erm i° n matrix 

^^fV^A (2-9) 

\ 2f) ls w i^uvy 

Using these discretizations, we have simulated the improved and unimproved (without surface 
terms) models applying a combination of Fourier accelerated HMC with higher-order integrators. 
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Figure 2: Left: Reduced improvement term AS/N for different lattice sizes: 9x9 (squares), 15 x 15 
(triangles) and 25 x 25 (circles). Colors depict A = 0.8 (red), 1.0 (green), 1.2 (blue), 1.5 (magenta). 
Right: MC history of improvement term and fermion determinant (SLAC improved, N = 1 5 x 1 5, mi att = 0.6, 
A = 1.4 (green), 1.7 (red), 1.9 (blue)). 



3. Limitations of the Nicolai improvement 



For simulations of the improved model 
including dynamical fermions the expec- 
tation value of the bosonic action is inde- 
pendently of A fixed to 

(5 B ) =N = # lattice points. (3.1) 

Nevertheless the improvement term AS = 
L.v \ W x (d<p) x + W x (d<p)xJ does not nec- 
essarily vanish. Therefore, we analyze 
AS with SLAC fermions at different cou- 
plings and for different lattice masses m\ att 
= m/N s (Fig. |, left panel). Here N s de- 
notes the number of lattice points in spa- 
tial direction. In the continuum limit (m\ atl 
— ► 0) the improvement term consistently 
vanishes for every A. 



19- 



Figure 3: Mode analysis of ensembles in the physical 
(green, A = 1 .4) and unphysical (red, A = 1 .7) phase. Here 
p is the distribution function for the modulus of the lattice 
momentum averaged over 25,000 configurations (SLAC im- 
proved, N = 15 x 15, rai att = 0.6). 
For (AS) /N > 14% the behavior of improvement term and fermion determinant changes sig- 
nificantly (Fig. ||, right panel). The improvement terms dominates the bosonic action by more than 
one order of magnitude while (Sb) = N is still preserved. Additionally the fermion determinant 
grows drastically and so hinders the system from returning into the original region of configuration 
space. This instability can be explained by reconsidering the improved action 



5 B = W 2(d<p) x + W x 2 = £ (2 (d$) x (d(p) x + \ \W X 



+ AS. 



(3.2) 



This action allows for two distinct behaviors of the fluctuating fields. The physically expected 
behavior consists of small fluctuations around the classical minima of the potential. Alternatively, 
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A = 0.2 
A = 0.4 



Figure 4: Masses for bosons (<pi, <f>2, statistics 10 6 -10 7 configs) and fermions (statistics 10 4 configs) for 
improved (left) and unimproved (right) model with standard Wilson fermions. 



(3.2) allows for large fluctuations of kinetic and potential term to be compensated by the improve- 
ment term of opposite sign. In this situation, it is definitely no longer possible to extract meaningful 
physics. 

Analyzing the distribution of the fields in momentum space at A = 1.4 and A = 1.7 (Fig. |f) 
shows that for too large couplings A (or lattice masses «ii att ) the simulation samples only unphysical 
UV dominated configurations. Therefore at strong coupling a careful analysis of the improvement 
term during the simulation must be ensured in order to achieve reasonable simulations. 



4. Weak coupling 

In the regime of weak couplings (A < 0.4) we are able to match bosonic and fermionic masses so 
as to observe how well supersymmetry effects (e.g. the degeneracy of masses) are realized on the 
lattice. Furthermore continuum extrapolations of the different discretizations are compared to the 
result of continuum perturbation theory at one-loop order. 

4.1 Signs of supersymmetry at finite lattice spacing 

In an unbroken supersymmetric theory bosonic and fermionic masses coincide. In the lattice formu- 
lation the supersymmetry is broken explicitly (at least partially). This induces a possible breaking 
of the mass multiplets which is explored at different lattice spacings for A G {0.2,0.4}, m = 15 
(Fig. Even with a statistics of up to 10 7 configurations the masses of bosons and fermions can 
not be distinguished in our simulations. Additionally improved and unimproved models give the 
same results (within error bars) for A < 0.4. 

This demonstrates that for Wilson type fermions in a region where the simulations do not 
show unphysical UV effects the Nicolai improvement is not necessary. A stable simulation with 
the unimproved model is likely to provide the same results, at least in the continuum limit. 

4.2 Continuum extrapolation 

For the free theory the lattice masses can be computed analytically. To make contact with perturba- 
tion theory which is carried out in the continuum it is crucial to get a stable continuum extrapolation 
even for the interacting case. Extrapolations from finite lattice spacing to the continuum using stan- 
dard and twisted Wilson fermions for the improved model (m = 15, A = 0.3) are shown in Fig. || 
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Figure 5: Left: The continuum extrapolation of fermionic masses for X = 0.3 for the improved Wilson and 
twisted Wilson model. Here, the SLAC result is given for one single lattice size. For comparison the exact 
results for the free theory are also shown. 

Right: Continuum extrapolation of fermionic masses for the weakly coupled regime in comparison to the 
perturbative result. 



(left panel). These are based on lattice sizes N s € {20, 24, 32, 48, 64} and demonstrate that both for- 
mulations yield the same continuum result. Additionally this result also coincides with a prediction 
by the SLAC model on a finite lattice (N = 45 x 45). Furthermore the twisted Wilson fermions are 
much closer to the continuum limit than standard Wilson fermions for finite lattice spacing. There- 
fore our analysis of the intermediate coupling case uses only the twisted type of Wilson fermions 
and SLAC fermions. 



4.3 Comparison to perturbation theory 

For small X we compare the perturba- 
tive one-loop result for the renormalized 
mass 



m 1 1 1 



4A 2 \ 
3^37 + 



(4.1) 



to the continuum extrapolation of the lat- 
tice data (Fig. ||, right panel). All dif- 
ferent formulations are seen to coincide 
with perturbation theory. Even for the 
unimproved model with Wilson fermions 
the correct (supersymmetric) continuum 
limit is reached within error bars. 

5. Intermediate coupling 
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Figure 6: Masses of the improved and unimproved model 
with SLAC fermions on a 45 x 45 lattice and continuum ex- 
trapolated results for twisted Wilson fermions are compared 
with the perturbative one-loop result in the continuum. 



To explore the limitations of the one-loop calculation we have performed simulations with X G 
[0, 1.2] (see Fig. |6|). The continuum extrapolations of Wilson type fermions are only applicable up 
to X < 0.7 using lattice sizes of N s < 64 due to the improvement problems. To cope with this, we 
instead use SLAC fermions which allow for a much larger X range on the accessible lattice sizes. 
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For A > 0.6 the improved and unimproved model with SLAC fermions give slightly different 
results on a 45 x 45 lattice. To check which model is closer to the continuum limit additional 
simulations with N = 63 x 63 at X = 0.8 have been performed (Tab. [I]). The data unveil that the 
unimproved model suffers from stronger finite a effects. Therefore the correct continuum limit is 
reached for both SLAC models but the improved SLAC model is closer to the continuum limit on a 
finite lattice. 

6. Conclusions and outlook 

We have performed a detailed analysis of the Nicolai im- 
provement in the Wess-Zumino model. This improvement 
introduces new problems due to the sampling of unphysi- 
cal (high-momentum) states. Additionally with high statis- 





improved 


unimproved 


45 


10.22(26) 


11.49(9) 


63 


10.54(15) 


10.70(19) 



tics bosonic and fermionic masses can not be distinguished Table 1: Fermionic masses for the 
for Wilson type fermions at finite lattice spacing in the weak SLAC derivative on two different lat- 

. . . j. . ,. c , • a a tice sizes for X — 0.8. 

to intermediate coupling region for both the improved and 

unimproved formulations. Even without improvement the correct continuum limit is reached. 
Therefore the term "improvement" is somewhat misleading. Only for SLAC fermions in the inter- 
mediate coupling region the improved action is closer to the continuum limit on finite lattices. 

More detailed results on this model including the discussion of discrete symmetries, the ab- 
sence of finite size effects and the effects of negative fermion determinants can be found in [|fl. 

In order to access the region of stronger couplings (X > 1.5) further algorithmic improve- 
ments are necessary. With the help of the PHMC or RHMC algorithm and improved solvers and 
preconditioners we are confident to obtain strong coupling results in the near future. The elaborate 
algorithms will then be used to explore the Jf = 1 Wess-Zumino model in two dimensions where 
a spontaneous supersymmetry breaking is expected and to study supersymmetric CP^ models in 
two dimensions. 
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